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Abstract — A newr DC-SQUID is proposed that exploits the 
properties of the grain boundary junctions in high Tc super- 
conductors. The orientations of the grain boundaries are 
chosen in such a way to estabhsh a tt/2 (rather than or 
7r) phase difference between the equilibrium phases of the 
two Josephson junctions in the SQUID loop. This property 
is exploited to increase the sensitivity and direction depen- 
dence of the SQUID for measuring magnetic flux without 
additional flux generating coils. 

Keywords — SQUID, High Tc superconductivity, Quasiclas- 
sical theory. 

I. Introduction 

IN a usual Josephson tunnel junction, the current / pass- 
ing through the junction is related to the phase difference 
(f between the two superconductors by m 



- Junction 



I — Ic aimp. 



(1) 



The equilibrium phase difference is zero in the absence of 
current (current is also zero when ip = n, but this corre- 
sponds to the maximum of the Josephson energy rather 
than a minimum). The simple relation (III) should be re- 
placed with a more complicated one 0, B in the case of 
constriction junctions, but the equilibrium phase difference 
will still he a,t if = 0. 

A conventional DC-SQUID consists of a superconduct- 
ing ring with two ordinary Josephson junctions on its op- 
posite arms. The total current passing through the SQUID 
is the sum of the currents crossing each of the Josephson 
junctions. As a result of interference between the two cur- 
rents, the total current-phase relation depends on the flux 
$ threading the ring. In particular, the critical current (the 
maximum current allowed before going to non-stationary 
state with finite voltage drop across the junction) is given 

by|] 

Ic == IcQ COS ( ^ ) , (2) 



' where (j)e is related to the external flux $ threading the 
loop through 0e = 27r$/$o, with <I>o = h/2e being the 
flux quantum, h the Planck constant, and e the charge of 

, electron. Here we neglect the self-inductance of the loop. 
Notice that (pe in Eq. (||) has period 2n which corresponds 
to one flux quantum. At small 0e, the relation (^ is al- 
most flat {Ic oc 0g); consequently, the sensitivity of the 
measurement of small fluxes is low. Moreover, the critical 
current does not depend on the direction of the magnetic 
field. To avoid this problem, one can apply a biasing flux 
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Fig. 1 

Symmetric and Asymmetric d-wave grain boundary junctions, with 

45° misorientation of the order parameters on both sides, are used 

as two Josephson junctions of the DC-SQUID. 



to move the working point of the SQUID away from the 
flat region. However, this biasing flux increases the noise 
as well as the complexity of the device and sometimes can 
have unwanted influence on the measured system. 

Therefore, it is desirable to shift the working point of 
the SQUID away from the flat region without applying 
an external magnetic field. This is possible by replacing 
one of the Josephson junctions with a junction with non- 
zero equilibrium phase difference. A DC-SQUID with a tt- 
phase difference between the equilibrium phases of the two 
Josephson junctions has already been studied Q. Here we 
propose another type of SQUID (we call it a 7r/2-SQUID) 
with the equilibrium phase difference of 7r/2 between the 
junctions. The frustration caused by the two junctions pro- 
duces the desired shift of the working point. In the next 
section we introduce such a device using high Tc super- 
conducting grain boundary junctions. In section III we 
study the 7r/2-SQUID using a simple model. Section IV 
is devoted to the quasiclassical calculation of the current 
in the d-wave structure and demonstrating agreement with 
the simple model of section II. Finally, we summarize our 
results in section V. 

II. High Tc 7r/2-SQUID 

In high Tc superconductors, because of the d-wave sym- 
metry of the order parameter, the current-phase relations 
are nontrivial and generally depend on the orientations of 
the order parameters on both sides of the junction Iq], 0. 
If we keep only the first two harmonics, we can write the 
current-phase relation as 



I = Ii sin if ~ I2 sin 2ip. 



(3) 



When I2 > /1/2 the equihbriuin phase wiU occur at 
if = ifo, where (po is neither nor tt. In particular, when the 
first harmonic vanishes, the equilibrium phase difference 
happens at ipo = 7r/2. These types of junctions, known 
as 7r/2-junctions, can be realized at the grain boundary 
junction between two d-wave superconductors with crys- 
talline orientations of 0° and 45° with respect to the grain 
boundary. At such a junction (known as an asymmetric 
junction), there is an exact cancellation of the first har- 
monic and the second harmonic dominates the current g , 
M. These junctions have been realized and tested experi- 
mentally m. A symmetric junction can also be made with 
45° misorientation angle by choosing -22.5° and 22.5° crys- 
talline orientations at the two sides. This junction is very 
different from the asymmetric one; specifically, the exact 
cancellation of the first harmonic does not occur in this case 
[|6|. In a realistic junction with non-ideal transparency and 
roughness at the boundary, the second harmonic gets sup- 
pressed by the boundary imperfections more strongly than 
the first harmonic and, as a result, the latter dominates the 
current. This leads to a conventional junction (O-junction) 
with the usual current-phase relation of Eq. (|). These 
junctions have also been studied experimentally 8|. 

By combining symmetric and asymmetric junctions, we 
can have a SQUID with 7r/2 equilibrium phase difference 
between the junctions. Fig. 1 illustrates such an structure. 
First we study this system using a simple model. The more 
realistic numerical simulation of the system using a quasi- 
classical model will come after that. 

III. Simple Model 

Consider a DC-SQUID, consisting of a conventional and 
a 7r/2 Josephson junction, with current-phase relations 

h = Id sin (fi, l2 = -Ic2 sin 2ip2 , (4) 

respectively. To simplify the calculations, we will assume 
Id — 2/c2 = 2/co. The more general case of arbitrary 
Ici/Ic2 can also be treated as the solutions are similar. We 
also neglect the self inductance of the ring in our calcula- 
tions. When there is a flux $ threading the SQUID ring, 
the phases at two junctions are related by fi — (p2 — <Pe- 
The total current passing through the SQUID is then given 

by 

/ = /co [2 sin((/j2 + </»£) -sin 2(^2]- (5) 

The critical current is the maximum current in this current- 
phase relation. To find the critical current, one has to find 
the solutions to the equation dl/d(p2 = 0. These solutions 
are (pe + 2-Kn and (— (/>e + 27rn)/3, with n being an integer 
number. Substituting into Eq. (|5|) we find the maximum 
current to be 



/. 



ico Max s sm 



{4>e + mi) 



(6) 



2 " 1^ 3^ 

Here "Max" means the maximum of the three cases with 
n = 0,±l. 

Fig. 2 shows the flux dependence of the critical current 
using Eq. (||). Notice that the periodicity of the func- 
tion is TT (half a flux quantum) although the period of the 
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Fig. 2 

Solid line is the critical current as a function of external 

flux for. the simplified model. dashed line shows the 

sine-functions in Eq. (6). 



argument of the sine-functions in Eq. (o) is actually Btt. 
The TT-periodicity of Eq. (o) is a result of the symmetry 
/(0e + ■"■, '■P2) = ^(0e: V2 + tt) in Eq. (P). In other words, 
changing the applied flux by half a flux quantum shifts the 
current-phase relation by tt and therefore does not change 
the critical current. This is evidently different from the 
ordinary SQUID for which the period is 27r (one flux quan- 
tum). The important feature here is that the 0e = does 
not coincide with the maximum of the function as it does 
in usual the DC-SQUIDs. Therefore, an increase in sensi- 
tivity for measuring small fluxes is gained. Moreover, the 
direction of the fiux is also detectable by measuring the 
critical current. 

IV. QUASICLASSICAL CALCULATIONS 

As was mentioned, the current-phase relations in d-wave 
grain boundary junctions does not usually have the simple 
forms of Eq. (Q) . One has to find the realistic current-phase 
relation using a microscopic theory, taking into account the 
imperfections of the boundary. Here we use a quasiclassical 
model to study this system. The equations we solve are the 
Eilenberger equations 



d ^ ^ ^ 

VF-irG'^(vF,r) + [cj?3 + A(vF,r),G'^(vi.,r)] =0, (7) 
or 

where 



A=( °. ^V G.(vF,r)=r^J l^^ ]. (8) 



At 



A is the superconducting order parameter and G^j (v^ , r) is 
the matrix Green's function, which depends on the electron 
velocity on the Fermi surface vj?, the coordinate r, and the 
Matsubara frequency to = (2n + l)7rT', with n being an 
integer number and T the temperature. We also need to 
satisfy the normalization condition 



1 - /-/. 



(9) 
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Fig. 3 
Current-phase relation obtained from the quasiclassical 

MODEL USING r = 1.4. 



Fig. 4 
Critical current versus the external flux for different 

VALUES OF r. 



In general, A depends on the direction of vp and is deter- 
mined by the self-consistency equation 



A(vf , r) = 27rN{0)T ^ {V{vp, Vp)Uv'p, r))^^ (10) 

LJ>0 

where V{vF,'v'p) is the interaction potential. Solution of 
matrix equation M) together with (nw determines the cur- 
rent density j(r) in the system 



j(r) 



-47rie7V(0)r^ 



(vFffa;(Vi=^,r)), 



(11) 



w>0 



In two dimensions, A^(0) = m/27r is the 2D density of states 
and (...) = J {d9/2TT)... is the averaging over directions of 



the 2D vector vp. 

The transparency of the boundary is described by the 
parameter Do |1^, ||ll|]. To incorporate the roughness, we 
assume a narrow scatterin g re gion of width d between the 
two superconductors [Q, [Q. The roughness is therefore 
parameterized by p = d/l, where I is the mean free path of 
the quasiparticles in this scattering region. In our calcula- 
tions, we use Dq — p — 0.5 which are reasonable values for 
a realistic system. Different widths, Wi and W2, are cho- 
sen for the two arms and the ratio of the widths is denoted 
by r = W2/W1. The details of our numerical method for 
solving these equations is given in ||ll|. 

Fig. 3 displays the current-phase relation at different val- 
ues of the external flux 0e for a SQUID structure with 
r — 1.4. At (j)e = 0, the graph has two local maxima, 
which is the signature of having a mixture of first and sec- 
ond harmonics. As 0e is increased, the two maxima become 
closer and at some value of (j)e, they will be equal. At this 
point, the global maximum jumps from one local maximum 
to another. This actually explains the sharp kinks at the 
minimum positions in Fig. 2 (and also in Fig. 4). 



Fig. 4 shows the critical current as a function of (f>e for 
different values of r. One can immediately see the simi- 
larities between these graphs and the one obtained from 
our simple model in Fig. 2. Specifically, the maxima again 
occur with period of tt, although the 7r-periodicity is not 
exact here (the real period is 27r). The maximum of the 
critical current is shifted from (jje — 0, as was desired. One 
can also see from Fig. 4 that the maximum range of vari- 
ation of Ic is obtained at r = 1.4. This is actually the 
optimal ratio of the widths for SQUID design. 

V. Conclusions 

We proposed and studied a DC-SQUID which has tt/2 
difference between the equilibrium phases at its two junc- 
tions. The 7r/2 phase difference is easily obtained using the 
properties of high Tc grain boundary symmetric and asym- 
metric junctions. The frustration caused by this phase dif- 
ference moves the working point of the SQUID to a point at 
which the maximum critical current is obtained at a non- 
zero external flux. This will increase the sensitivity of the 
SQUID at small magnetic field, with a dependence on the 
direction of the flux, without any external biasing coil. 

We introduced a simple model and compared its results 
with the ones obtained from a microscopic quasiclassical 
model. The agreement between the two models is ac- 
ceptable. We found that the 7r/2-SQUID exhibits other 
new features, different from a conventional SQUID or a 
TT-SQUID, among those is the $o/2 (instead of $0) peri- 
odicity of the critical current. We find the optimal ratio 
between the widths of the two arms to be r = 1.4. 
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